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We show that when a suitable entanglement generating unitary operator depending on a parameter 
is applied on N qubits in parallel, a precision of order 2~ N in estimating the parameter may be 
achieved. This exponentially improves the precision achievable in classical and in quantum non- 
entangling strategies. 
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The Mandelstam-Tamm version of the time energy un- 
certainty relation [| and its rigorous developments Q 
form the basis of quantum enhanced methods for pa- 
rameter estimation such as those used in evolving fre- 
quency standards |3(- Giovannetti, Lloyd and Maccone 
|4| have recently done beautiful work to introduce a gen- 
eral framework to classify methods of such quantum en- 
hancement. A unitary transformation depending on the 
parameter to be estimated is applied to a suitably pre- 
pared system of N probes and then an appropriate ob- 
servable is measured. A separable initial state is labelled 
C and an entangled state Q. Similarly, measured observ- 
ables which are direct products of individual qubit ob- 
servables are designated C and those which introduce 
entanglement are designated Q. Thus the strategies of 
"parallel" measurement on N qubits are classified as CC, 
CQ, QC, QQ. Their result is that quantum metrology 
in the strategies CC and CQ only achieve the classical- 
precision limit of order e/ \/A , where e is the dispersion 
of measurement results for each probe; but the quantum 
metrology strategies QC and QQ can achieve a preci- 
sion of order e/N if the probes are prepared in a suitable 
entangled state. In contrast, for A sequential measure- 
ments on a single probe to achieve the same precision of 
order 1 /N requires a running time and hence duration of 
quantum coherence to be N times longer. 

We propose here a parallel strategy which exploits the 
Hilbert space of N probes more fully than in previous 
work and thereby attains an exponentially enhanced 
precision. In discussing parallel strategies previous work 
have only considered applying on the prepared probes a 
unitary operator which is a direct product of A unitary 
operators each acting on a single probe. We will show 
here that if instead we consider applying on the probes, 
an entanglement generating unitary operator U = e~ lBH 
which cannot be written as a direct product of one-probe 
operators, then we can obtain an exponentially enhanced 
precision in estimating the parameter 9. The fundamen- 
tal reason for this improvement is that there are an ex- 
ponentially large number of mutually commuting observ- 
ables for the iV-probe system whereas the number of mu- 
tually commuting single-probe operators is only of order 
N. For instance, for the A-qubit system each qubit can- 



not have more than one commuting observable and hence 
there are exactly A commuting observables of the form 
l®---®!®^®!®...! each of which acts non- 
trivially on only the jth qubit, with j = 1, 2, . . . , N. We 
construct sets of 2 N_1 mutually commuting Hermitian 
operators each of the form A 1 ® A 2 eg) • • • ® A N which acts 
non-trivially on all the A qubits. The unitary entangle- 
ment generating operator is chosen to be a product (but 
not a direct product) of the exponentials of these Hermi- 
tian operators times — i6. In this way, we fully exploit the 
quantum parallelism which is at the heart of exponential 
violations of local realism 

5] 

and of the well-known expo- 
nential speed-up achieved in certain quantum computa- 
tion tasks . For the probe system of A qubits we obtain 
the best possible precision of order e 2~ N dictated by the 
uncertainty principle for such entanglement-generating 
unitary operators. The previous precision limit of order 
e/N is the best possible for unitary operators which are 
direct products. We thus obtain a new characterization 
of eight parallel quantum metrology strategies as XYZ 
where each of X, Y and Z can be Q or C. See Fig.Q] Here 
X and Z specify presence or absence of entanglement in 
the probes and observables respectively as in and 
the new label Y specifies whether the unitary operator 
U applied on the probes is an entanglement-generating 
operator or a direct product. 

In our parallel strategy, a precision of order e2~ N is 
obtained using an entanglement generating operator that 
is a product of 2 JV_1 mutually commuting unitary oper- 
ators, but entanglement at the probe preparation stage 
and measurement stage are inessential. To obtain the 
same precision in the sequential strategy of Ref. one 
needs an exponentially large number of applications of 
a unitary operator on a single qubit; hence the time for 
which quantum coherence has to be maintained will be 
longer by a factor of order 2^. The increase in precision 
we achieve with respect to the parallel strategy of Ref. is 
however at the cost of a "Hamiltonian" or "generator" of 
U which might be more difficult to implement experimen- 
tally. Here again, quantum mechanics just might come to 
the rescue. We propose a quantum-optics model of laser 
light interacting with an A-qubit system, say a poly- 
atomic molecule, via a generalized Jaynes-Cummings in- 
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FIG. 1: The eight parallel strategies XYZ for estimation 
of a parameter occurring in an unitary operator U applied 
to N probes classified according to types of state prepara- 
tion,unitary operation, and measurement. Each of X, Y, Z 
takes values C or Q. Y — Q if U is an entanglement gen- 
erating operator (as considered in this paper) and Y = C 
otherwise. X, Z specify presence of entanglement (Q) or ab- 
sence of entanglement (C) in the prepared probe state and in 
the measured operator respectively. 



teraction which, in principle, could achieve the exponen- 
tially enhanced precision. The practical implementabil- 
ity of this specific model is at the moment unclear, but 
would be an interesting subject for further theoretical 
and experimental investigation. 

Quantum limits on precision: We recall first the 
Mandelstam-Tamm uncertainty relations. Suppose we 
apply an unitary operator U — e~ l9H to a probe state 
\4>) where H is a self-adjoint operator and 8 a parameter 
to be estimated. We obtain the state 



(i) 



and then measure the observable X on this state with a 
view to estimating the parameter 8. Schwarz inequality 
yields 



AX AH > \(ij(e)\[H ) X}/(2i)\m)\ 



no 



(2) 



where AX, AH denote the dispersions in X, H respec- 
tively. The resulting precision in estimating 8 is thus 
given by the uncertainty relation 



= AX/ 
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2AH 



(3) 



In the case when the parameter 8 is estimated using an 
estimator 8 0St and the estimation is repeated v times, the 
Cramer-Rao bound Q was used in Ref. to prove the 
generalized uncertainty relation for the error estimate 
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where (• • • ) av denotes statistical average. The quantum 
limits on precision are obtained by noting that 



AH < l(A n 



•^min) ; 



(5) 



where A max , A m ; n denote respectively the maximum and 
minimum eigenvalues of H which we assume to be finite. 
Hence, the Mandelstam-Tamm pj and Braunstein-Caves 
quantum limits on precision are given respectively by 



> 



(Amax A m in) 



and 



S v 6 > 



AmirOv^ 



(6) 



(7) 



It is clear that quantum parallel strategies to improve 
precision should aim to maximize AH on the iV-probe 
quantum state. 

Parallel strategies for N qubits: Consider first, as in 
|4( the operator H to be a direct sum of the operators 
ffl acting on the j-th probe, each ff J having the same 
maximum dispersion AH } < d. Then 



A0ff< Nd 



(8) 



i=i 



As noted in Ref. this maximum dispersion is in fact 
reached when we choose the iV-qubit state to be an 
equally weighted superposition of the eigenvectors of H 
with maximum and minimum eigenvalues. With A max = 
NXm and Amin = NX m where Am and A m are respec- 
tively the maximum and minimum eigenvalues of each 
H J we get the above equation with d — (Am — A m )/2. 
Note for purposes of comparison, that the dispersion of 
a sum of K classical variables H 3 ci with a factorized joint 
probability distribution is given by 



K 



\ 



E( A ^) 2 



(9) 



The maximum quantum dispersion is thus v N times 
the classical value when K = N and AH^ = d. This 
has been exploited in Ref. 0. We now show that the 
dispersion of H for N qubits can be made exponentially 
larger by appropriate choice of H. Consider, the operator 
identity 



JV 



= H + iA 



(10) 



where a 3 x and a 3 y are Pauli matrices for the j-th qubit, and 
H and iA denote respectively the Hermitian and anti- 
Hermitian parts of the operator on the left-hand side of 
the equation. Explicit expressions for H and A are con- 
veniently stated in terms of the matrices er J (±l) defined 
by 



^(+1) 



■1) 



(11) 
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We obtain 



H 



^2 H(n,r 2 ,...,r N ) 



N_ ^even 



where 



A' 



H(r 1 ,r 2 ,...,r N ) = (-l) N -/ 2 (gl^(r j ) 



and 



r 1) r 3 ,...,r JV =±l 
AT_=odd 



where 



and finally 



N 1 



(12) 



(13) 



A= A(n,r 2 ,...,r N ) , (14) 



A(r 1 ,r 2 ,...,r N ) = (-l)( N - 1 V 2 (g)ai(r j ) , (15) 



(16) 



in A. In contrast with the parallel strategy in Ref. [j, we 
now have 



AH < 2 



N-l 



AA<2 



N-l 



(21) 



Interestingly, and in contrast with Ref. the maximum 
dispersions of H and A are now reached in the separable 
iV-qubit states 

ITT •••!!), IU-..U), (22) 

each of which has 

AH = 2 iV " 1 , AA = 2 N - 1 , (23) 

where we have denoted the eigenstates of a z with eigen- 
values + 1 and —1 by | f ) and |J.) respectively. 

For comparison, the maximum classical dispersion for 
the sum of 2 7V ~ 1 classical variables with factorized prob- 
ability distribution, each variable having maximum dis- 
persion 1, would be 2( Ar ~ 1 )/ 2 . 

Parallel strategy CQC with exponentially en- 
hanced precision: We start from one of the factorized 
./V-qubit states given above and apply one of the unitary 
operators U — e~ l6H or U = e~ l6A given above. In par- 
ticular, we obtain 



is just the number of j's with Tj = — 1 or the number of 
<r 3 yS in the iV-fold product of Pauli matrices in H and 
A. Both H and A are sums of 2 N ~ 1 products of Pauli 
matrices, each product having eigenvalues ±1 and hence 
maximum dispersion 



N 



A(gK(r,)<l. 



(17) 



The standard anti-commutation rules between Pauli 
matrices leads to 



^ Oj )& 3 {r'i ) = (t j {r' } )a J (rj )r^ 



(18) 



for r,-,rj = ±1. Note that rir 2 ...r N = (-1) . The 
anti-commutation rules then imply that the set of 2 N ~ 1 
products of Pauli matrices occurring in H (or A) consti- 
tutes a set of mutually commuting observables. Hence 



-i8H 



n 



-i9H(r 1 ,r2,...,r N ) 



r 1 : r 2 ,.-.,r N = ±l 
JV_=even 



and 



iOA 



n 



^-i9A(ri,r 2 ,...,r N ) 



(19) 



(20) 



r\ ,r 2 ,-..,r N = ±l 
N_ —odd 



Further any of these 2 N 1 observables in H anti- 
commutes with any of the 2 A '~ 1 observables occurring 



\i> H {6)) = e- WH ITT---TT) 



(24) 



= 008(2^-^)1 T...T)-isin(2^-^)||... |); 



\M9)) 



-i0A 



| TT — TT> 



(25) 



cos(2 JV -^)|T...T)+sm(2 JV -^)U..4) 



We may now measure the probability that all qubits are 
in the up-state given by the expectation value of the di- 
rect product of projection operators 



X 



N 



1 



(26) 



We obtain, for example for the state \ipii(d)) 



(X) = i[l + cos(02*)] , AX = i|sin(02*)| . (27) 
Hence, the quantum precision of estimating 9 is given by 



= 2 



-N 



(28) 



which achieves the best allowed by the uncertainty re- 
lation since AH = 2 N ~ 1 . A similar precision is ob- 
tained by using the state 1^(0)). It is clear that we 
only needed to use the parallel strategy CQC (no entan- 
glement in the probe-preparation or measurement stage, 
but an entanglement-generating unitary operator applied 
on the probes) to obtain this exponential enhancement 
of precision. 
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A generalized Jaynes-Cummings model: We pro- 
pose a quantum optics model of interaction of laser light 
with N qubits which could, in principle, be used to obtain 
exponentially enhanced precision. We apply a unitary 
operator to the infinite-dimensional vector space which 
implies an entanglement generating 'operation' on the 
reduced density operator for N qubits. 

Suppose we shine laser light of frequency f2 « Ntv on 
an iV-atom molecule each atom of which can be approx- 
imated to be a two-level atom or qubit with level sepa- 
ration hw . We assume that on irradiation, the molecule 
does not dissociate but gets excited to a higher level in 
which each atom is excited to the higher level. The sec- 
ond quantized interaction Hamiltonian is assumed to be 



H = H + H! 



Ho 



N 

IjJ \ -> 

2 ^ 



a{ + ttata , (29) 



N 



h.c. 



i=i 



(30) 



Here h.c. denotes Hermitian conjugate, a and denote 
annihilation and creation operators for photons of fre- 
quency SI, \{o x + i<J y y denotes the level raising opera- 
tors for the j-th qubit and g is a real coupling constant. 
It is an exactly solvable Hamiltonian with eigenvalues 



1 Or 

2 7 2 



A± = Sl(n - 
where n = 0, 1, 2, . . . and 

= ^{ujN - SI) 2 + g 2 (n + 1)2 2N . 
The corresponding eigenstates are, 

i(ar±|n}|to}+Afc|n + l}|&>) > 

|to) = |T...T>, |<M = II •••!), 



(31) 



(32) 



(33) 



(34) 



a± = ^1±{luN -SV>/Sl 1 



(35) 



and a'c 



Interesting physics about Rabi oscil- 



lations between up and down qubit states may be read 
off from these equations even off-resonance, i.e., f2 ^ Nu>. 
For the present we specialize to SI = Nm which implies 
[H , Hi] = 0. We deduce that for r = and r = 1 



e- itH In 



_ -it(n+l/2)Nu 



x[cos(2 N - 1 d)\n + r)\(f> r ) 

-isin(2 Ar - 1 6»)|n + l-r)|0i_ r )] 



(36) 



where 9 — tgyn + 1. Thus the reduced density opera- 
tor for the N qubits obtained by tracing over the pho- 
ton states undergoes the positivity and trace preserving 
transformation 

p{t)= h r p(Q)K , h r = e- im \<j> r ){cj> r \, (37) 

r=0,l 

where H is the iV-qubit spin operator given by Eq. I|1(J|I . 
For measurement of the observable X of Eq. I|26() with 
(X) = tr p(t)X we obtain exactly the same results on 
precision as in Eqs. I|27l) and l|28l) of the previous sec- 
tion. Investigation of practical implcmcntability of this 
model for achieving exponentially enhanced precision is 
an important remaining task. 

In conclusion, we have shown that the parameter esti- 
mation associated with suitable entanglement-generating 
unitary operators may lead to an exponential enhance- 
ment of accuracy over both classical schemes and non- 
entangling quantum schemes. The exciting thing is that 
unlike quantum computation, quantum metrology might 
not face formidable problems of fighting decoherence if a 
suitable interaction Hamiltonian can be found. We have 
shown that the suggested strategy may be implemented 
in principle using a generalized Jaynes-Cummings inter- 
action. 
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